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We study how the stability of the fractional quantum Hall effect (FQHE) is influenced by the geometry of 
band structure in lattice Chem insulators. We consider the Hofstadter model, which converges to continuum 
Landau levels in the limit of small flux per plaquette. This gives us a degree of analytic control not possible in 
generic lattice models, and we are able to obtain analytic expressions for the relevant geometric criteria. These 
may be differentiated by whether they converge exponentially or polynomially to the continuum limit. We 
demonstrate that the latter criteria have a dominant effect on the physics of interacting particles in Hofstadter 
bands in this low flux density regime. In particular, we show that the many-body gap depends monotonically on 
a band-geometric criterion related to the trace of the Fubini-Study metric. 
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I. INTRODUCTION 
A. Overview 

The theoretical proposal and numerical observation of time- 
reversal breaking lattice models exhibiting fractional quantum 
Hall (FQH) phases, reviewed in Refs. 1 and 2, have sparked 
intense interest in these models (“fractional Chern insulators” 
or FCIs) as a promising venue for observing FQH physics 
without the need for large external magnetic fields. The ma¬ 
jority of our theoretical understanding of the the fractional 
quantum Hall effect (FQHE) has been framed in the context of 
continuum Landau levels, which occupy a special and highly 
atypical point in the space of all single-particle Hamiltonians. 
It is a matter of both theoretical and experimental interest to 
identify which of the special properties of Landau levels (LLs) 
are essential to the existence of the LQHE: as LCIs rapidly 
move towards experimental reality, one would like to easily 
identify experimental parameters where LQH-like phases are 
most robust. 

Relevant to these aims are several recent studies inves¬ 
tigating the role of the algebra of band-projected density 
operators in an LCI in shaping the excitation spectrum of 
LQH phases [3-7]. Intriguingly, these results point to a role 
for band geometry, i.e., non-topological properties of Chern 
bands, in determining the stability of LQH states in LCIs; in 
particular, the Fubini-Study metric describing the embedding 
of the occupied bands in the Hilbert space of all bands appears 
naturally when one considers a long-wavelength expansion of 
this operator algebra [7]. That reference gave qualitative ar¬ 
guments linking properties of this algebra with the stability of 
FQH-like states which were numerically investigated in [8]. 
The results obtained there suffered from the limitation that un¬ 
derlying band geometry could only be varied indirectly; more¬ 
over, one would like to be able to unify these numerical ex¬ 
amples with a theoretical picture capable of making analytic 
statements. 
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For these reasons, in this article we investigate the condi¬ 
tions proposed in [7] in the context of the Hofstadter model 
[9]. This provides an ideal laboratory for studying the inter¬ 
play between band geometry and FQH physics, due to the fact 
that it offers a controlled limit (that of vanishing flux per pla¬ 
quette) which converges to continuum LLs, which offer opti¬ 
mal conditions for FQH states (in a sense to be made precise 
below) and allows us to obtain analytic results not possible 
for generic FCIs. We note that, although the Hofstadter model 
involves a net magnetic flux, numerical studies [10, 11] have 
shown its FQH states may be adiabatically connected to those 
in FCIs with no net flux. A significant literature on the FQHE 
in the Hofstadter model exists [12-19]. The Hofstadter model 
is also a subject of ongoing experimental work; Hostadter 
bands have recently been realized with cold atoms in optical 
lattices [20-22] and in graphene superlattices [23, 24]. 

In what follows, we show that the band-geometric cri¬ 
teria proposed in [7] meaningfully quantify the stability of 
LQH phases as one perturbs the band geometry of the Hof¬ 
stadter model away from the continuum limit. We show 
that momentum-space fluctuations of the Berry curvature and 
Lubini-Study metric have exponential convergence to their 
LL limit, while a third, nontrivial condition related to the 
averaged trace of the Fubini-Study metric has slower, poly¬ 
nomial convergence. Numerical exact diagonalization stud¬ 
ies of several different FQH states show that the many-body 
gap depends monotonically on this third criterion, even when 
momentum-space fluctuations are negligible. This constitutes 
strong evidence that the band geometry plays a key role in 
the stability of the FQHE in settings where lattice effects are 
important. 

B. Geometry of Chern bands 

In this section and the next, we present the background ma¬ 
terial necessary to formulate our hypothesis for the stability 
of LQH-like states in the Hofstadter model and fix notational 
conventions which will be used in subsequent calculations. 

In momentum space, we write the band Hamiltonian of a 
lattice model as i7(k) = i?Q,(k)|k, a)(k, a|, where 
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a = 1,..., N indexes the bands. A band is parameterized by 
its Bloch function rt“(k), an A^-component vector in terms of 
which |k, a) = b), where |k, b) is the Fourier 

transform of the &th basis orbital. Topological order in a band 
is indicated by nonzero values of a corresponding topologi¬ 
cal invariant, which for two-dimensional systems with broken 
time-reversal symmetry is the first Chem number 


Cl = ABz{Ba) I^TT, (1) 

where Abz is the area of the Brillouin zone (BZ), (• • •) de¬ 
notes the average over the BZ, normalized so (1) = 1, and 
the Berry curvature [25, 26] of the band a in terms of Bloch 
functions is 
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One may additionally define a “quantum metric” on the BZ by 
introducing the Fubini-Study metric on the band Hilbert space 
[27] and using (k) to pull it back to the BZ, yielding 
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We emphasize that, despite their mathematical origin, the cur¬ 
vature and metric are quite natural physical quantities: observ¬ 
able effects of Berry curvature are reviewed in [28], and the 
quantum metric is experimentally accessible through current 
noise measurements [29]. 

In the present work, we consider three geometric condi¬ 
tions for the stability of FCIs involving Ba and g°‘ which were 
first proposed in Ref. 7 (in a refinement of previous arguments 
[1, 5]). These are that (i) the Berry curvature (2) is constant 
over the BZ, (ii) all components of the quantum metric (3) are 
constant over the BZ, and (iii) the following quantity 

i9(k) = det<7“(k)-^^^^^|^>0 (4) 

vanishes everywhere, i.e. the metric determinant inequality 
(4) is saturated. In Ref. 7 a similar inequality involving the 
trace of the metric 


T(k)=trp“(k)-|B„(k)|>0 (5) 

was proved, and it was shown that T(k) = 0 automatically 
implies 79 (k) = 0. 

Conditions (i), (ii) and (iii) were not chosen arbitrarily: the 
main result of Ref. 7 is that they are sufficient to guarantee that 
the set of band-projected density operators in an FCI is closed 


under commutation, and obeys a generalization of the Girvin- 
MacDonald-Platzman algebra [30] arising in the continuum 
FQHE. Density-density interactions are essential to the stabil¬ 
ity of any FQH phase, so the more accurately the algebra is 
reproduced (specifically, as measured by the conditions), the 
more stable FQH-like phases should be: we refer to this as the 
“band geometry hypothesis.” 

The validity of the band geometry hypothesis has been 
tested for a variety of FQH-like phases appearing in several 
different FCI models in [8]. We take a different approach in 
the present work: we consider a sequence of lattice models 
(the Hofstadter model at </) = 1/N, defined below) which con¬ 
verges to continuum FQHE physics as N ^ oo, which lets us 
quantitatively investigate how the stability of EQH states de¬ 
pends on band geometry in a way that was not possible in 
Ref. 8. 


C. Hofstadter model in the limit of small flux per plaquette 


The quantum dynamics of an electron on a two-dimensional 
square lattice in the presence of a transverse magnetic field 
was first studied by Harper [31], and its self-similar band 
structure was further elucidated by Azbel [32]. With these 
caveats, for brevity we refer to the following tight-binding 
Hamiltonian as the Hofstadter model [9]: 


H = -t 


(ij) 


clc,exp(2„^ 


A ■ d£] + h.c. 


( 6 ) 


with i,j indexing sites on a square lattice with lattice constant 
a. In the Landau gauge, A = Bx y and all eigenstates have 
trivial y dependence ip{x = na,y) = with the re¬ 

maining X dependence obeying the discrete Harper equation 

Ipn+l + 1pn-l + 2 C0s(27r(/)n - ky)i>n = (7) 

where e is the energy in units of t and (j) is the magnetic flux 
per lattice cell in units of the flux quantum do = hje. For 
rational d = P/Q, (7) is invariant under n ^ n + Q and one 
may apply Bloch’s theorem over the enlarged magnetic unit 
cell consisting of Q x 1 lattice plaquettes. 

In [19], the Hofstadter model was studied perturbatively for 
d near rational values. We recall some results relevant to the 
case of interest here, d = 1/A^: the W —>■ oo limit of the 
Hofstadter model is the continuum Landau level (LL) problem 
in a precise sense. We first recall that in the Landau gauge, 
states in the LL problem are spanned by two commuting free 
boson algebras generated by a,a^,6, where the first pair 
act as ladder operators for the LLs and the second pair raise 
and lower momentum within a LL. A major result of [19] is 
that Hofstadter bands may be written as superpositions of LLs 
by a unitary change of basis, given for large N by 
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D. WKB Hofstadter wavefunctions 


In this section we recall details of the Hofstadter eigenfunc¬ 
tions in the WKB approximation [19, 33, 34] (not all of which 
have been published) which are needed for our calculation in 
Sec. IIB. In the large-limit, we may regard the lattice site 
index n in the Harper equation (7) as a continuous variable 
X = n/N. Replacing discrete differences by derivatives, we 
obtain a second-order differential equation in which permits a 
WKB approximation: 

-I- [(e -I- 2) -I- 2 cos(27ra;)] -i/Tfa:) = 0. (9) 

Here we have set ky = 0, since increasing ky is equivalent to 
a translation in x in the Landau gauge. 

The classical turning points of (9) occur at x = mzL f] for 
integer m and /3 = (I/Stt) cos“^ (—e/2 — 1); in the large- 
limit /3 ^ 1/ \/2TrN, so that in the large-W limit the two first- 
order turning points atm — /3, m + ^ coincide at the minimum 
of the potential to form a single second-order turning point (so 
that the region corresponding to the oscillatory WKB solution 
vanishes). We use a quadratic approximation to the cosine 
potential near the turning points and use the WKB solution 
only in the classically forbidden region. 

We define the Bloch wavefunction piecewise in three re¬ 
gions of the magnetic unit cell, with 


AU{a,0 + BV{a,0, 


uu{x) = < 


0 < X < a/3; 

C-iptxpix) + D-ijj-pix), 

afi < X < 1 — a/3; 

F[/(a,^-e(l))+Gna,e-e(l)), 

1 — a/3 < X < 1. 


The first and third regions are in the vicinity of the second 
order turning point (integer x), with the solution given by 
parabolic cylinder functions U{a,^), V{a,^) as discussed in 
Ref 19, where a = ^ (e-f 4) and /'(x) = 2^/ttNx are 
suitably rescaled variables. The second region is the classi¬ 
cally forbidden region, where the solution decays exponen¬ 
tially. There is no canonical point at which the wavefunctions 
from different regions should be matched, which is expressed 
by a dimensionless free parameter a > 1. Its value does not 
change the qualitative nature of our conclusions, and we use 
the value a = 2.3 adopted in Ref. 19 so that our results are 
compatible with that work. 

The Bloch function in the second, classically forbidden re¬ 


gion is [34] 


Uk,ii(a;) = 


\/sinh^(x) 


— ik^xN^ 
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( 10 ) 


u+{x-ky/2 t:) 

(correcting a misprint in Ref. 19), where u±{x) is 


u± (x) = exp 


±N [ p{y)dy 
Jll2 


( 11 ) 


and for future convenience we have defined 


p{x) = cosh ^ [—e/2 — cos(27rx)]. (12) 


We have the symmetries u_(x) = m+( 1 — x) and m_{x) = 
l/u_|_(x). Effects which are nonperturbative in A^ in (10) and 
what follows are governed by powers of e~°'^, where we de¬ 
fine 


rl-P 

tf = / p{x)dx 

Jp 
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with the asymptotic value given by 4/7r times Catalan’s con¬ 
stant, G = "'^31 follows we only 

keep contributions of 

Finally, in (10) A is the normalization of the entire Bloch 
function u^. This is dominated by contributions from the first 
and third regions, where jukP has its maxima [34]: 


'^W\/TErf(a)' 

Contributions from region II enter \A\‘^ at 0{N~'^), but we’ve 
found that the coefficient is numerically negligible at a = 2.3, 
and we drop it in the rest of the analysis to avoid the compli¬ 
cations involved in attempting to satisfy the self-consistent set 
of equations that would result. 

Note that in (10) the ky dependence has been re-inserted 
“by hand”, under the assumption of continuous translational 
dependence. As a consequence, we cannot recover any infor¬ 
mation about fluctuations in the ky direction and must appeal 
to symmetry arguments presented in section IIA to constrain 
the full functional dependence from the kx behavior. 


E. Summary 

Our results are grouped into two sections. In Section II, we 
And analytical expressions for aspects of the quantum metric 
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of the lowest Hofstadter band in the limit of small flux per 
plaquette (p. Considering lattice effects as perturbations of the 
lowest Landau level in this limit, we calculate the average of 
the metric over the Brillouin zone and the k dependence of 
fluctuations about this average. We also provide nonperturba- 
tive results, in particular statements about the k dependence 
of the metric from symmetry considerations. The results we 
obtain in Section II together with similar results for the Berry 
curvature found in Ref 19 allow us to describe the behavior 
of the determinant and trace inequalities. 

In Section III, we study multiple FQH-like states by nu¬ 
merically diagonalizing many-particle interactions projected 
to the lowest Hofstadter band and finding the corresponding 
many-body spectra. Our main findings are that geometric 
quantities can be distinguished by their functional dependence 
on N in the continuum limit, vanishing either exponentially 
or polynomially, and that the latter of these appear to have the 
dominant effect on the stability of fractional phases as mea¬ 
sured by the size of the many-body gap. Our conclusions are 
given in Sec. IV. In Appendix A, we present explicit calcula¬ 
tions for N = 3 which illustrate the results of Sec. II. 


II. QUANTUM GEOMETRY OF THE HOFSTADTER 
MODEL 

A. Symmetries of the quantum metric 

In this section we prove results [Eqs. (23) and (24)] con¬ 
cerning the symmetry of the quantum metric of the Hofstadter 
model over the Brillouin zone. This lemma is needed in the 
following section, where we use symmetry properties to in¬ 
fer the full k dependence of the quantum metric from the 
dependence of the WKB wavefunctions. The argument be¬ 
low closely parallels that given for the Berry curvature in Ap¬ 
pendix B of Ref. 19, and all of the results are readily verifiable 
numerically. 

We proceed by comparing eigenfunctions for the Hofstadter 
Hamiltonian written in two different Landau gauges, = 
Bxy and = —BySc. As indicated, in this section we 

identify quantities computed in the two different gauges by the 
superscripts (Y) and (A), respectively. The Harper equation 
in the (V) gauge is 

+ [2 cos(27r(/)a; - ky) - e] = 0, (15) 

where In ths (A) gauge, the 

Harper equation takes the form 

'(/’if+i + + [2 cos{2-K(j>y + k^) - e] ^ = 0, (16) 

where instead the x dependence is trivial: %p^^\x^n'a) = 
. Equation (15) is transformed to (16) by the clock¬ 
wise rotation {x, y) —>■ (—y, x), {kx, ky) —>■ {—ky, kx), so this 
also sends eigenfunctions to eigenfunctions. Since the corre¬ 
sponding Bloch functions also transform trivially, this means 
that the behavior of the curvature and metric is determined 


entirely by the transformation properties of the momentum 


derivatives involved in their definition: 

9Z\kx,ky) = g^y^\-ky,kx), (17) 

9xy\kx, ky) = —gxy'^ {—ky, kx), (18) 

9m\kx,ky) = gi^\-ky,kx), (19) 

B^^\kx,ky)=B^^\-ky,kx). (20) 


<Y) 

We observe that if ipP is a solution of (15) with k = ky, 
then (with the index n interpreted cyclically modulo N) 
is also a solution with the same value of e and k = ky + 2TT(f>. 
Taken together with the periodicity imposed by the magnetic 
unit cell, this implies that B^^'^ and each component of g^^^ 
has periodicity 27r/A in both kx and ky, for any (j> of the form 
M/N. A parallel argument holds for the (A) gauge. 

In order to compare the gauges (A) and (A), we consider 
a real-space unit cell of N x N sites, so that each eigen¬ 
state in both gauges appears with an artificial A-fold degen¬ 
eracy. Ordinary lattice translations by N sites in the x and y 
directions commute with each other and either Hamiltonian, 
so we may define Va Ak) as the A-dimensional subspace 
spanned by eigenstates corresponding to the band a, with 
momentum eigenvalues k now defined in the reduced BZ of 
size 2tt/N x 2tt/N. We may obtain a basis for Va Ak) for 
fixed k by taking the eigenvectors {kx + 27rm/A, ky) for 

m = 0,1,..., A — 1. 

The curvature and components of the metric for the sub- 
space Va ' (k) may be defined as the sum of the correspond¬ 
ing quantity over all elements of an orthonormal basis for that 
space, and it readily follows from the definitions (2), (3) that 
this is independent of the choice of basis (in the literature, this 
is sometimes referred to as a “gauge invariance” of the Berry 
curvature and quantum metric, which should not be confused 
with the electromagnetic gauge symmetry under discussion 
here.) Using the Bloch functions corresponding to the basis 
vectors %p^^\kx + 21:771/q, ky), we have 

g,y[V^^\k)]=Ng^’^^\kx,ky), (21) 

B[V^^\k)]=NB2'\kx,ky). (22) 

Corresponding statements can be made for the (A) gauge. 
Now let Uyx be the unitary operator mapping the Hamilto¬ 
nian in the (A) gauge to the Hamiltonian in the {Y) gauge. 
This sends the subspace Va^\k) to lAl^^(k), so we may use 
(k)) as a basis for (k) to And gij (k)] = 
gij\y^'^{k)\, and likewise for the curvature. Together with 
(21), (22) and the 2t: / N x2t: / N periodicity of all quantities in 
question, this proves the transformation properties (17)-(20) 
apply within the gauge (A), i.e. with (A) on the right-hand 
sides of (17)-(20) replaced by (A). In particular, 

iYg{kx,ky) =tig{-ky,kx), (23) 

deig{kx,ky) = deig{-ky,kx)- (24) 

In the following section, we will use a WKB analysis in 
the (A) gauge to show that the leading-order kx dependence 
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of B, Qxx and gyy is cos{Nkx)- The symmetry (23), (24) 
then constrains the full momentum dependence of B, tr g 
and det p to be a function of cos{Nkx) + cos{Nky) and 
cos{Nkx) cos{Nky). One may proceed to show that, of these 
two forms, the curvature depends only on the sum [19]. Oddly, 
despite having the same symmetry as B, we have found that 
tr g and det g depend on both the sum and product of cosines. 
This may be seen from numerics, or established analytically 
for small N using the results of Appendix A. 


B. Quantum metric fluctuations 

As explained above, the band geometry hypothesis predicts 
that the two most important factors influencing the stability of 
the FQHE in a lattice system are fluctuations in the Berry cur¬ 
vature and quantum metric, as a function of k. The N depen¬ 
dence of the amplitude of these fluctuations may be computed 
by using the WKB approximation to the Bloch wavefunction 
of the lowest Hofstadter band in the definitions of B and g 
[Eqs. (2) and (3)]. This was done for curvature fluctuations in 
Ref 19, but the results for fluctuations of the quantum metric 
we present below are new. 

Given the Bloch function (10), we can take derivatives and 
calculate the curvature and metric using Eqs. (2) and (3), with 
the inner product over the band index replaced by an integral 
over the unit cell and N f dx. To lowest order, the 

parabolic cylinder function solution corresponds to continuum 
Landau levels, for which the metric and curvature are uniform. 
We therefore make the assumption that the leading-order os¬ 
cillatory contribution to the band geometry comes solely from 
the WKB piece of the wavefunction in the second region, 
which we described in detail in Sec. 1D. 

Contributions from the second term of (3) contain two fac¬ 
tors of and are neglected in our analysis. Differentia¬ 

tion of u±{x — ky/2'K) with respect to ky brings down factors 
of Np{x — ky/2'x)/2'K, while the terms with oscillatory be¬ 
havior in kx that we seek arise only from differentiation of 
the explicit exponential factors in (10) with respect to kx- As 
a consequence, if we choose to measure the fluctuations of 
a quantity X{kx,ky) via X{kx) = X{kx,0) — Ar(0, 0), all 
terms involving \u±{x)\'^ cancel at We obtain 

~gxx = \A\’^N^e-^^h cos{Nkx), (25) 

gxy = 0 + Oi\A\'^N^e-^^^), (26) 

gyy = \A\^N^e-^^h cos(Wfc,), (27) 

B = 2\A\‘^N^e-^^l2 cos{Nkx), (28) 

where A is given by (14), a is given by (13) and the needed 
integrals are 


p{xY dx 
ap (27r)2 sinhp(x) ’ 


(31) 


where we have used the fact that all integrands are even under 
reflection about x = 1/2 to fold the region of integration. The 
integrand obtained for gxy is, instead, odd about x = 1/2 and 
vanishes; the only contributions to gxy come from terms we 
neglect at our level of approximation. 

In Refs. 19 and 34, [sinhp(x)]“^/^ was approximated by 
its average value (which is 0(1)) and factored out of all 
integrations. This is a rather severe approximation, since 
[sinhp(x)]“^/^ diverges as x —)■ /3, where relevant N de¬ 
pendence enters; in fact, under the change of variables u = 
sin^(7rx) we And 


r.l/2 


h = 2 


dx 


sinhp(x) 




^ (ttch/S) 


du-^ log 
du 


.1 — \/\ — V? 


where the lower limit of integration goes to 0 as /2N. The 
integrand diverges as u —>■ 0, and we have Jq ~ (l/27r) log N. 
Note that this divergence is irrelevant in practice, as it is dom¬ 
inated by 6“*^^ wherever it appears. 

The other integrands are progressively better behaved as 
X ^ jd, and their integrals are all finite as N ^ oo. Nu¬ 
merical integration at a = 2.3 gives 


h 

h 

h 


0.114- 

0.125- 

0.0214 


0.251 

0.068 

(32) 

Vn ^ 

N ’ 

0.292 

0.063 

(33) 

Vn ^ 

N ’ 

0.002 


(34) 

N 



We And excellent agreement with the amplitude of the 
cos{Nkx) + cos{Nky) dependence, measured numerically, of 
each quantity (red points in Eig. 1). Unlike the Berry curva¬ 
ture, the metric components have additional momentum de¬ 
pendence which is not of this form (see Sec. II A), but we 
And that the asymptotic estimates obtained in this discussion 
also accurately describe the total RMS fluctuation of gxx and 
gyy (black points in Eig. 1), despite the fact that, due to the 
choice of Landau gauge, our analysis is unable to reconstruct 
the full momentum dependence of the components of the met¬ 
ric for large N. In the case of gxy, we simply plot the leading 
N dependence \A\'^N'^e~'^‘^^ as a guide to the eye, since we 
cannot determine the constant of proportionality at the cur¬ 
rent level of approximation. Similar remarks hold for fluctu¬ 
ations in H; we And its N dependence is well described by 

(l/20)|A|4w5e-‘^^. 


h = 2 



x(l — x) 
sinhp(x) 


dx, 


(29) 


C. BZ-averaged metric components 


p(x) dx 
a/3 (27r) sinhp(x)’ 


In this section, we provide details on the perturbative cal¬ 
culation of BZ-averaged band-geometric quantities. The final 
condition appearing in the band geometry hypothesis involves 
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FIG. 1. (color online). Magnitude of Brillouin zone fluctuations 
of band-geometric quantities, as a function of N. Red points cor¬ 
respond to the amplitude of the cos{Nk:,;) + cos{Nky) term in the 
Fourier series, while black points correspond to the RMS fluctuation 
of each quantity found from numerical integration over the BZ. Blue 
lines show the approximation given in Eqs. (25)-(28). Our analysis 
only predicts the form of the leading N dependence for g^y and D, 
which is shown with a thin line (the overall proportionality constant 
being undetermined.) 


the Brillouin zone average of the determinant (or trace) of the 
quantum metric, as opposed to the fluctuations calculated in 
the previous section. Since this does not require knowing the 
k dependence of these quantities, they may be computed more 
simply by changing to the Landau level basis using (8), rather 
than by manipulating the WKB wavefunctions directly. 

Using (8), the Hofstadter ground state |0) = U^|0), where 
|0) is the lowest Landau level, is given in the LL basis \m) 
(m > 0) as 

1 TT^ 1 TT^ 2419 TT^ 

9-217]^ ■ 

/- / 1 TT 1 TT^ 371 TT^ \ i,\ 

+ 12V™(^^ + ,..)|8) + ..,, (35) 

where we have only written out the terms needed for comput¬ 
ing matrix elements to Note that we are able to work 

to fourth order despite having only computed the exponent of 
W to third order, because the term with the fewest boson op¬ 
erators appearing in the exponent of W is always — a'^): 




this is a manifestation of the rotational symmetry of the con¬ 
tinuum being broken by the square lattice. It can be verified 
that adding fourth-order corrections to the coefficients in (8) 
do not affect any of the terms in (35), or any of the results 
below. 

We now seek to compute BZ averages of metric compo¬ 
nents, using 


= Re tr [P^n (1 - Pg) fyPg]. (36) 

where Pg = |0)(0| projects onto the lowest band of the Hof¬ 
stadter model, which may be expressed in terms of LL pro¬ 
jectors using (35). Though the form of the operator W (and 
hence |C))) was derived in the Landau gauge, the perturba¬ 
tive corrections we incorporate preserve the degeneracy of the 
Landau levels and therefore are independent of the choice of 
gauge at our level of approximation. This point is discussed 
in more detail in Sec. VI A of Ref. 19. The calculation of the 
averaged metric is easiest to carry out in the symmetric gauge, 
where the quantum-mechanical position operators may be ex¬ 
pressed in terms of the LL ladder operators a, and angular 
momentum raising and lowering operators bfbas 

X = ^{a — — ib — ib^), (37) 

y = ^{a + + ib — ib^). (38) 


In (36) these operators appear sandwiched between the 
ground-state projector and its complement. This simplifies 
the problem considerably: because b, b^ do not mix LLs, no 
term containing them contributes in (36), and we may neglect 
the angular momentum degree of freedom and work purely 
in the space of LL indices (as the notation of (35) suggests). 
Furthermore, because |()) to 0{N~^) only has support on LLs 
|0), |4) and |8), terms of the form aa and vanish as well. 
The only matrix element that contributes to (36) is 


~ 1 TT^ 1 TT^ 

^ ^ 647^ 


it. 


311 


TT 


9 • 211 ^4 


-FO(At-5) ^ 


(39) 


Eq. (39) suffices to compute the averages of all elements 
of gij as (gy;;^) = {gyy) = 1/2 + (0|a'l'a|0) and (g^y) = 
0 + 0{1/N^). Although gxx and gyy do not have the same k 
dependence, symmetry requires that their BZ averages be the 
same. 

Because the BZ average of the curvature is always exactly 
quantized to the Chern number and all higher moments of all 
quantities (considered as distributions over the BZ) vanish to 
any order in 1 /N (due to the fact that our analysis of the WKB 
wavefunctions in Sec. IIB showed that the amplitude of all k 
dependence falls off exponentially in N), it follows that 


(T) = 2 {gxx) - 1 

1 TT^ ^ 1 TT^ ^ 311 TT"! 

(D) = {gx.f - 1/4 

1 TT^ ^ 1 TT^ ^ 313 TT^ 

"" 96i\^'^647^'^9-2iii^'^"" 


( 41 ) 
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Unlike the momentum-space fluctuations of the metric 
components and Berry curvature, the trace and determinant 
inequalities have a much slower asymptotic decay which is 
polynomial in 1/A^ (Fig. 2), rather than exponentially small 
in N. Therefore, at large N these conditions should be the 
dominant factors dictating the stability of FQH-like phases. 
In practice, this regime occurs for N larger than ^ 15, as can 
be seen in Fig. 2. 



FIG. 2. (color online). Plot of the BZ-averaged value of the trace (5) 
and determinant inequalities (4). Points are obtained from numeri¬ 
cally integrated values of (T) and (D). The blue lines are given by 
Eqs. (40), (41). The thin black lines show the fluctuation amplitudes 
T, D plotted in the last row of Fig. 1 for comparison. 

Finally, as an overall check on the consistency of our large- 
N approximation, one may eliminate (gxx) from (40), (41) to 
obtain 4(11) — (T)((T) -F 2) = 0; the left-hand side of this 
relation is plotted in Fig. 3. All neglected terms, namely (gxy) 
and BZ fluctuations of terms quadratic in B and g, vanish ex¬ 
ponentially with the N dependence expected from the WKB 
calculation of Sec. IIB. 


III. STABILITY OF FQH STATES 
A. Scaling of many-body gaps 

We investigate the degree to which the above band- 
geometric criteria influence the stability of an FQH state by 
adding repulsive interactions and numerically computing the 
many-body gap for different values of (j> = 1 /N. Lowering 



FIG. 3. (color online). Plot of the residuals in the implicit relation¬ 
ship between (40) and (41). Points are obtained from values of (T) 
and (D) numerically integrated over the BZ. The blue line is a plot 
of the predicted N dependence, , as a guide to the eye. 


the flux per plaquette changes the relative strength of the in¬ 
teraction, since the size of the system (size of the magnetic 
unit cell) increases while the number of particles remains con¬ 
stant. In order to compare gaps at different N, the strength of 
the interaction potential must be scaled with N as follows. 

In the spirit of the WKB analysis, let n = 1,..., W be 
a wavefunction defined on the lattice and V'(a^) be its large- 
continuum approximation, with x = n/N. Under the (non¬ 
trivial) assumption that ipn has support over almost all tight- 
binding orbitals, the requirement that both ipn and tpi^) be 
normalized to unity dictates the scaling tpn '/N’tp{x). Con¬ 
sideration of the matrix elements of a two-body (5-function in¬ 
teraction V{x) = Vq5{x) shows that 

(■0^, U |^/>^, = Vo J dxtp^{x)*'ijj^{x)*'ij}^{x)ijj^{x) 

n 

so that the strength of the corresponding discrete on-site in¬ 
teraction should scale as Vdisc VqN. Similar considerations 
show that for the two-body nearest-neighbor repulsion used to 
stabilize the fermionic Laughlin state and for the three-body 
delta-function interaction stabilizing the Moore-Read state, 
the leading scaling should be Vdisc ~ VqN'^. Because the 
single-particle bandwidth vanishes exponentially fast with N, 
the many-body gap should be scaled in the same way, and 
we define Ajc = NA for the bosonic Laughlin state and 
Asc = A^A for the fermionic Laughlin and Moore-Read 
states, where A is the gap obtained with a fixed interaction 
strength. Similar scaling arguments can be used to compare 
FCI models with different numbers of bands [8]. 


B. Many-body computations 

In order to study the influence of band geometry on FCI 
phases, we carried out exact numerical diagonalization of a 
repulsive interaction Hamiltonian projected to the lowest Hof- 
stadter band for Np = 8 particles partially filling the lowest 
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band. To study the bosonic Laughlin state at z/ = 1/2, we used 
a lattice of 4 x 4 unit cells and a two-body on-site repulsion; 
for the fermionic Laughlin state at z/ = 1/3 we used a 4 x 6 lat¬ 
tice with a two-body nearest-neighbor repulsion, and for the 
bosonic Moore-Read state at = 1 we used a 2 x 4 lattice 
with a three-body on-site repulsion. For geometric reasons, 
we used system sizes of the form N = rn? up to = 169 for 
the bosonic Laughlin state, N = Qm? up to iV = 294 for the 
fermionic Laughlin state, and N = 2m? up to iV = 162 for 
the bosonic Moore-Read state. 

For each system, we verihed that the many-body ground 
state had nontrivial topological order corresponding to the ap¬ 
propriate FQH state via properties of energy and entanglement 
spectra. We required that the many-body spectrum exhibit a 
quasidegenerate ground state, with the gap A to excited states 
larger than the spread in ground-state energies. We compute 
the particle entanglement spectrum by tracing out four par¬ 
ticles from the density matrix formed by an equal superpo¬ 
sition of all ground states, and required that this spectrum 
be gapped, with the counting of eigenvalues below the gap 
in each momentum sector given by the appropriate counting 
rules [35, 36] for the corresponding FQH state. 

In Fig. 4 we plot the scaled gaps as a function of (T), the 
Brillouin zone average of the trace inequality (5). We note a 
clear correlation between an increasing value of the gap and 
increasing saturation of the inequality, with the gap continu¬ 
ing to increase even for values of N for which BZ fluctua¬ 
tions of B and g are negligible (as shown in Fig. 1). Apart 
from a change in horizontal scale, dependence of the scaled 
gaps on {D) is visually identical; comparing (40) with (41) 
shows that, in the large-7V limit, {D) « 2{T) up to terms of 
0(iV-4). 


IV. DISCUSSION 

The Hofstadter model provides an ideal laboratory for 
studying FCI phenomena; the existence of a controlled limit 
in which its spectrum converges to continuum LLs, which 
makes the model amenable to perturbative expansion in 1 /N 
and allows a controlled study of the relationship between band 
geometry and FQHE physics. In this paper, we have shown 
that there is a natural distinction between effects that are non- 
perturbatively small in 1/7V, such as the BZ fluctuations of the 
Berry curvature and quantum metric, and effects which have 
a perturbative expansion in 1/iV, such as the trace and deter¬ 
minant conditions. The behavior of the many-body gaps ob¬ 
tained via exact diagonalization indicates that the latter effects 
dominate for large N. In the future, it would be interesting to 
study models where the role of the determinant inequality can 
be isolated. 

The Hofstadter model has been realized experimentally in 
optical lattices of cold atoms [20-22] and graphene superlat¬ 
tices [23, 24]. The geometrical criteria for the stability of 
Chern bands are readily computable single-particle properties 
that act as a meaningful proxy for the many-body gap, and 
hence may serve as a useful guide in selecting experimental 
parameters and couplings which are most favorable to the ex- 
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FIG. 4. (color online). Scaled many-body gaps as a function of BZ- 
averaged trace inequality (5) for FQH-like states at various values of 
N. Insets show behavior of points at larger values of N. Error bars 
reflect the finite accuracy of our diagonalization code. (a). Laughlin 
state of Np = 8 bosons atv— 1/2. (b). Laughlin state of Np = 8 
fermions at z/ = 1/3. (c). Moore-Read state of Np = 8 bosons at 
z/ = 1. 


istence of a stable FQH state. In this sense, the present work 
complements Ref. 8, in which the effects of band geometry on 
the many-body gap were measured by varying the couplings 
of several FCI models. In that study, BZ fluctuations of the 
Berry curvature were found to have the largest influence on 
the gap, with the trace condition a subdominant factor. 

In this paper we have presented a perturbative treatment of 


































9 


the band geometry of the Hofstadter model at </) = 1/7V in the 
limit N —> oo, but our analysis may be extended to the most 
general case considered in [19], namely (p = P/Q + M/N 
with N ^ M. Our own single- and many-particle computa¬ 
tions suggest that the case p = M/N is qualitatively similar 
to the results discussed above, with all quantities taking values 
which interpolate among those obtained at M = 1. Different 
physics is encountered when perturbing around nonzero P/Q, 
since the lowest band now has Chern number ci = Q > 1. 
FCI states in bands with ci > 1 may be mapped onto multi¬ 
layer FQH states [37], albeit with significant distinctions [38- 
40]. The specific case of the Hofstadter model (and general¬ 
izations) at (/) = P/Q has been examined in [18, 41-43]. 

One could also extend our analysis to excited bands of the 
Hofstadter model, corresponding to higher LLs. Finally, ob¬ 


taining an analytic formula for the quantitative dependence of 
the gap (or other many-body observables) on band-geometric 
parameters remains an important open problem. 
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Appendix A: Closed-form band geometry for N = 3 

In [44], an analytic expression was obtained for the Berry 
curvature of an arbitrary three-band Hamiltonian. In this sec¬ 
tion, we obtain analogous expressions for components of the 
quantum metric, for the purpose of determining what state¬ 
ments may be made about the full momentum dependence of 
the quantum metric of the Hofstadter model at general N. 

We follow the notational conventions of Ref. 44. Just as 
an arbitrary two-band Hamiltonian may be written in terms 
of Pauli matrices and its eigenstates described using a Bloch 
spin, we may write any three-band Hamiltonian as i?(k) = 
a(k) -f b(k) • A, with the projector onto the ath band given 
by Pa = |(1 + y/Srialk) ■ A), for a = 1,2,3. Here A 
is an eight-component vector, the entries of which are the 
Gell-Mann matrices forming the fundamental representation 
of SU{3). These obey 

{^abc ^fabc'j^c- (Al) 

Using these structure constants we may dehne antisymmetric 
and symmetric vector products on eight-component vectors 
x,y as (x X y)a = fabc^bYc and (x * y)a = V^dabc^byc, 
respectively. 

Barnett, Boyd and Galitski parameterize the relationship 
between H and Pa via Ua = ‘CQ(7ab -I- b * b), with 

= 72 -V|2 ’ 7 a =2|b|cOS0„, 

where 


where we take a = \ for the lowest band. The fact that this 
expression involves the dot product instead of the antisymmet¬ 
ric product means that convenient cancellations arising in the 
computation of the Berry curvature do not happen here. We 
find 


12[2cos(20) + l]4 

X |5 -I- 16cos^ [1 -I- 5cos(20) -f cos(40)] 

— 8-\/2 cos 0 [1 -I- 2 cos(20)]^ cos{3kx) 

-I- [5 -I- 4cos(20)] cos( 6 fca;)} , 

with 6 = 6*1 (k) given by (A2). The expression for gyy is 
analogous, with the explicit factors of appearing in (A 6 ) 
replaced by ky, so that 


2 

“ 3[2cos(26») -f 1]4 

X {6 -I- 4cos(20) — cos(40) 

— [4cos(20) -I- 5] cos(3fca;) cos(3fcy)} . 


(A7) 


The momentum dependence of tr^ involves both the sum 
cos{Nkx) + cos{Nky) (through the k-dependent variable 9) 
and the product cos{Nkx) cos{Nky). We now show that the 
same situation holds for det g: using 


9xy — 


[4cos(20) -I- 5] sin(3A:a;) sin(3fcy) 
6[2cos(26l) -f 1]4 


(A 8 ) 


6'a = 2 cos ^ 


b • b * b 

IbP 


27r 


(A2) 


and, for clarity, we’ve suppressed k dependence of all quan¬ 
tities. The fact that Pa is a projector implies |nQ,p = 1 and 
riQ * ria = Hq,; the fact that it projects onto an eigenstate of 
H means [Pa,H] = 0, which implies b x Hq, = 0. 

The b vector corresponding to the N = 3 Hofstadter 
Hamiltonian in the Landau gauge is given (up to an irrele¬ 
vant permutation of rows and columns) following Eq. (4) of 
Ref 44: 

b(k) = |cosA: 2 ,,sinA:a;, ^ sinfcj^ — | cos fcy, cos fcx, 

— sin kx , cos kx , sin kx , ^ sin ky + ^ cos ky |, 

which is of constant length |bp = 6, implying b • db = 0. 
Additionally, 

b • b * b = 3^'^^[cos(3fc3;) -f cos(3fcy)], (A3) 


therefore, 6a and the other coefficients only depend on mo¬ 
mentum through the function cos(3fca;) -f cos{3ky). This is 
the cos{Nkx) + cos{Nky) momentum dependence found for 
the Berry curvature in Ref. 19. 

Components of the quantum metric may be computed as 

^ tr {Padk,Po^dk^ Po) (A4) 


we have 


detp = 


1 


36[2 cos(26») -f 1]6 
X {29 - 12cos(26i) - 48cos(46») 

— 12cos(60) — 2 cos( 8 d) 

- 1-6 [2 cos(20) -f 1] cos(3fca;) cos(3fcy)} . 


(A9) 


As proved in Sec. II A, this is the most general form of k 
dependence a function with the lattice symmetries of trp 
and det g is allowed to have. Despite being fully consis¬ 
tent with that proof, this result was somewhat unexpected: 
the Berry curvature has the same lattice symmetries, is de- 
hned as a similar combination of k derivatives of Bloch func¬ 
tions, and enters into the band geometry on an equal footing 
as the quantum metric [7], yet only depends only on the sum 
cos{Nkx) -b cos{Nky). 

In Fig. 5 we plot the curvature and each component of the 
metric for (j) = 1/3, where the fluctuations are largest in mag¬ 
nitude. The lattice symmetry properties proved above are evi¬ 
dent. 


(A5) 
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FIG. 5. (color online). Plots of the metric components and Berry 
curvature over a minimal periodic cell of the BZ (of size 2n/N x 
2tv/N) for the Hofstadter model at A'’ = 3. 



